(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     98992,       2594]
NotebookOptionsPosition[     91156,       2456]
NotebookOutlinePosition[     91536,       2473]
CellTagsIndexPosition[     91493,       2470]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\:7b2c5\:7ae0\:7b2c\:4e8c\:8bd5\:9a8c\:573a", "Title",
 CellChangeTimes->{{3.784971016355488*^9, 
  3.784971026684877*^9}},ExpressionUUID->"8961cf75-fdc5-4109-b21b-\
7a5dddeaa7ee"],

Cell[BoxData[
 RowBox[{"Clear", "[", "GuassJordon", "]"}]], "Input",
 CellChangeTimes->{{3.785016057965209*^9, 3.7850160625309753`*^9}},
 CellLabel->"In[49]:=",ExpressionUUID->"499f085a-2c33-40be-ba39-97a7fef4f643"],

Cell[BoxData[
 RowBox[{
  RowBox[{"GuassJordon", "[", "B_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "A", ",", "AA", ",", "K", ",", "F", ",", "a", ",", "S", ",", "b", ",", 
      "c", ",", "x", ",", "y", ",", "k", ",", "i", ",", "j", ",", "m", ",", 
      "n", ",", "p"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"A", "=", "B"}], ";", "\[IndentingNewLine]", 
     RowBox[{"c", "=", 
      RowBox[{"First", "[", 
       RowBox[{"Dimensions", "[", "A", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"F", "=", 
      RowBox[{"IdentityMatrix", "[", "c", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"i", "=", "1"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"i", "<", 
        RowBox[{"c", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"a", "=", 
           RowBox[{"Max", "[", 
            RowBox[{"Abs", "[", 
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"Table", "[", 
                 RowBox[{"j", ",", 
                  RowBox[{"{", 
                   RowBox[{"j", ",", "i", ",", "c"}], "}"}]}], "]"}], ",", 
                "i"}], "]"}], "]"}], "]"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"b", "=", 
           RowBox[{
            RowBox[{"First", "[", 
             RowBox[{
              RowBox[{"Position", "[", 
               RowBox[{
                RowBox[{"Abs", "[", 
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{
                    RowBox[{"Table", "[", 
                    RowBox[{"j", ",", 
                    RowBox[{"{", 
                    RowBox[{"j", ",", "i", ",", "c"}], "}"}]}], "]"}], ",", 
                    "i"}], "]"}], "]"}], "]"}], ",", "a"}], "]"}], "//", 
              "Flatten"}], "]"}], "+", 
            RowBox[{"(", 
             RowBox[{"i", "-", "1"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}]}], "}"}], "=", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}]}], "}"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}]}], "}"}], "=", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}]}], "}"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"n", "=", 
           RowBox[{"i", "+", "1"}]}], ";", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"A", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "All"}], "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "All"}], "]"}], "]"}], "/", 
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "i"}], "]"}], "]"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"F", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "All"}], "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"F", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "All"}], "]"}], "]"}], "/", 
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "i"}], "]"}], "]"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"S", "=", "F"}], ";", "\[IndentingNewLine]", 
          RowBox[{"While", "[", 
           RowBox[{
            RowBox[{"n", "<", 
             RowBox[{"c", "+", "1"}]}], ",", "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"k", "=", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"n", ",", "i"}], "]"}], "]"}], "/", 
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "i"}], "]"}], "]"}]}], ")"}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{
                RowBox[{"A", "[", 
                 RowBox[{"[", 
                  RowBox[{"n", ",", "All"}], "]"}], "]"}], "=", 
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"n", ",", "All"}], "]"}], "]"}], "-", 
                 RowBox[{"k", "*", 
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "All"}], "]"}], "]"}]}]}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{
                RowBox[{"F", "[", 
                 RowBox[{"[", 
                  RowBox[{"n", ",", "All"}], "]"}], "]"}], "=", 
                RowBox[{
                 RowBox[{"F", "[", 
                  RowBox[{"[", 
                   RowBox[{"n", ",", "All"}], "]"}], "]"}], "-", 
                 RowBox[{"k", "*", 
                  RowBox[{"F", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "All"}], "]"}], "]"}]}]}]}]}], "}"}], 
             ";", 
             RowBox[{"n", "++"}]}]}], "]"}]}], "}"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"i", "++"}]}]}], "\[IndentingNewLine]", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"K", "=", "A"}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", 
      RowBox[{"c", "-", "1"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"AA", "=", "A"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"p", ">", "0"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"A", "[", 
            RowBox[{"[", "p", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"A", "[", 
             RowBox[{"[", "p", "]"}], "]"}], "-", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"A", "[", 
                RowBox[{"[", "q", "]"}], "]"}], "*", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"AA", "[", 
                  RowBox[{"[", 
                   RowBox[{"p", ",", "q"}], "]"}], "]"}], "/", 
                 RowBox[{"AA", "[", 
                  RowBox[{"[", 
                   RowBox[{"q", ",", "q"}], "]"}], "]"}]}], ")"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"q", ",", 
                RowBox[{"p", "+", "1"}], ",", "c"}], "}"}]}], "]"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"F", "[", 
            RowBox[{"[", "p", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"F", "[", 
             RowBox[{"[", "p", "]"}], "]"}], "-", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"F", "[", 
                RowBox[{"[", "q", "]"}], "]"}], "*", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"AA", "[", 
                  RowBox[{"[", 
                   RowBox[{"p", ",", "q"}], "]"}], "]"}], "/", 
                 RowBox[{"AA", "[", 
                  RowBox[{"[", 
                   RowBox[{"q", ",", "q"}], "]"}], "]"}]}], ")"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"q", ",", 
                RowBox[{"p", "+", "1"}], ",", "c"}], "}"}]}], "]"}]}]}]}], 
         "}"}], ";", "\[IndentingNewLine]", 
        RowBox[{"p", "--"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Return", "[", "F", "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{3.785027719708577*^9},
 CellLabel->"In[71]:=",ExpressionUUID->"39468ca8-ce8d-4ece-893c-13e29e110352"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Inverse", "[", "B", "]"}], "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.785027735496086*^9, 3.7850277440781183`*^9}},
 CellLabel->"In[73]:=",ExpressionUUID->"2139ec7d-9357-4f15-8c94-601722cbade6"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "1"}], 
      FractionBox["1", "3"], 
      FractionBox["1", "3"]},
     {"0", 
      RowBox[{"-", 
       FractionBox["1", "6"]}], 
      FractionBox["1", "3"]},
     {"2", "0", 
      RowBox[{"-", "1"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.785027744549819*^9},
 CellLabel->
  "Out[73]//MatrixForm=",ExpressionUUID->"b1e4a0d4-affe-4465-a7ab-\
1fba7e15a44e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"GuassJordon", "[", "B", "]"}], "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.7850160785017643`*^9, 3.7850160866333375`*^9}},
 CellLabel->"In[72]:=",ExpressionUUID->"125fc4cd-f5fb-4769-b181-a7e7a26ef470"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", 
       FractionBox["1", "2"]}], "1", 
      FractionBox["1", "4"]},
     {"0", 
      RowBox[{"-", "2"}], "1"},
     {"1", "2", 
      RowBox[{"-", 
       FractionBox["3", "2"]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.7850160870946665`*^9, 3.785027726595169*^9},
 CellLabel->
  "Out[72]//MatrixForm=",ExpressionUUID->"82492a5e-7ab0-49a8-934a-\
53fbd03ef124"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"GuassJordon", "[", "B_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "A", ",", "AA", ",", "K", ",", "F", ",", "a", ",", "S", ",", "b", ",", 
      "c", ",", "x", ",", "y", ",", "k", ",", "i", ",", "j", ",", "m", ",", 
      "n", ",", "p"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"A", "=", "B"}], ";", "\[IndentingNewLine]", 
     RowBox[{"c", "=", 
      RowBox[{"First", "[", 
       RowBox[{"Dimensions", "[", "A", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"F", "=", 
      RowBox[{"IdentityMatrix", "[", "c", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"i", "=", "1"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"i", "<", 
        RowBox[{"c", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"a", "=", 
           RowBox[{"Max", "[", 
            RowBox[{"Abs", "[", 
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"Table", "[", 
                 RowBox[{"j", ",", 
                  RowBox[{"{", 
                   RowBox[{"j", ",", "i", ",", "c"}], "}"}]}], "]"}], ",", 
                "i"}], "]"}], "]"}], "]"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"b", "=", 
           RowBox[{
            RowBox[{"First", "[", 
             RowBox[{
              RowBox[{"Position", "[", 
               RowBox[{
                RowBox[{"Abs", "[", 
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{
                    RowBox[{"Table", "[", 
                    RowBox[{"j", ",", 
                    RowBox[{"{", 
                    RowBox[{"j", ",", "i", ",", "c"}], "}"}]}], "]"}], ",", 
                    "i"}], "]"}], "]"}], "]"}], ",", "a"}], "]"}], "//", 
              "Flatten"}], "]"}], "+", 
            RowBox[{"(", 
             RowBox[{"i", "-", "1"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}]}], "}"}], "=", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"A", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}]}], "}"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}]}], "}"}], "=", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"b", ",", "All"}], "]"}], "]"}], ",", 
             RowBox[{"F", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "All"}], "]"}], "]"}]}], "}"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"n", "=", 
           RowBox[{"i", "+", "1"}]}], ";", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"A", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "All"}], "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "All"}], "]"}], "]"}], "/", 
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "i"}], "]"}], "]"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"F", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "All"}], "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"F", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "All"}], "]"}], "]"}], "/", 
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "i"}], "]"}], "]"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"S", "=", "F"}], ";", "\[IndentingNewLine]", 
          RowBox[{"While", "[", 
           RowBox[{
            RowBox[{"n", "<", 
             RowBox[{"c", "+", "1"}]}], ",", "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{
                RowBox[{"A", "[", 
                 RowBox[{"[", 
                  RowBox[{"n", ",", "All"}], "]"}], "]"}], "=", 
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"n", ",", "All"}], "]"}], "]"}], "-", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{"n", ",", "i"}], "]"}], "]"}], "/", 
                    RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{"i", ",", "i"}], "]"}], "]"}]}], ")"}], "*", 
                  RowBox[{"A", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "All"}], "]"}], "]"}]}]}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{
                RowBox[{"F", "[", 
                 RowBox[{"[", 
                  RowBox[{"n", ",", "All"}], "]"}], "]"}], "=", 
                RowBox[{
                 RowBox[{"F", "[", 
                  RowBox[{"[", 
                   RowBox[{"n", ",", "All"}], "]"}], "]"}], "-", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{"n", ",", "i"}], "]"}], "]"}], "/", 
                    RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{"i", ",", "i"}], "]"}], "]"}]}], ")"}], "*", 
                  RowBox[{"F", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", "All"}], "]"}], "]"}]}]}]}]}], "}"}], 
             ";", 
             RowBox[{"n", "++"}]}]}], "]"}]}], "}"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"i", "++"}]}]}], "\[IndentingNewLine]", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"K", "=", "A"}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", 
      RowBox[{"c", "-", "1"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"AA", "=", "A"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"p", ">", "0"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"F", "[", 
            RowBox[{"[", "p", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"F", "[", 
             RowBox[{"[", "p", "]"}], "]"}], "-", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"F", "[", 
                RowBox[{"[", "q", "]"}], "]"}], "*", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"p", ",", "q"}], "]"}], "]"}], "/", 
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"q", ",", "q"}], "]"}], "]"}]}], ")"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"q", ",", 
                RowBox[{"p", "+", "1"}], ",", "c"}], "}"}]}], "]"}]}]}], ";", 
          RowBox[{
           RowBox[{"A", "[", 
            RowBox[{"[", "p", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"A", "[", 
             RowBox[{"[", "p", "]"}], "]"}], "-", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"A", "[", 
                RowBox[{"[", "q", "]"}], "]"}], "*", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"p", ",", "q"}], "]"}], "]"}], "/", 
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"q", ",", "q"}], "]"}], "]"}]}], ")"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"q", ",", 
                RowBox[{"p", "+", "1"}], ",", "c"}], "}"}]}], "]"}]}]}]}], 
         "}"}], ";", "\[IndentingNewLine]", 
        RowBox[{"p", "--"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Return", "[", "A", "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.784971132509079*^9, 3.7849711815774403`*^9}, {
   3.784971270764412*^9, 3.7849713256287184`*^9}, {3.784971477866038*^9, 
   3.7849714783876452`*^9}, {3.784971581956236*^9, 3.7849717132237186`*^9}, {
   3.7849718351290417`*^9, 3.784971967979576*^9}, {3.7849721039424906`*^9, 
   3.7849721116329327`*^9}, {3.7849721791888475`*^9, 3.784972302579177*^9}, {
   3.7849723341743603`*^9, 3.7849724029945374`*^9}, {3.784972508707*^9, 
   3.784972509255528*^9}, {3.7849726371719847`*^9, 3.784972653206111*^9}, {
   3.784972683815285*^9, 3.7849727148320494`*^9}, {3.784972826717886*^9, 
   3.7849728272085743`*^9}, {3.7849729520538845`*^9, 
   3.7849729742973833`*^9}, {3.784973063594164*^9, 3.784973081978061*^9}, {
   3.78497325140307*^9, 3.7849732519007053`*^9}, {3.7849732872297716`*^9, 
   3.7849732877873197`*^9}, {3.7849733202388616`*^9, 3.784973321624193*^9}, {
   3.7849734092197285`*^9, 3.784973409436113*^9}, {3.7849734571628747`*^9, 
   3.7849734581991043`*^9}, {3.7849735965827417`*^9, 3.784973626021047*^9}, {
   3.7849737603284855`*^9, 3.7849738000851717`*^9}, {3.7849739806530075`*^9, 
   3.784973988069229*^9}, {3.7849740718292785`*^9, 3.7849741344550514`*^9}, {
   3.784974219932563*^9, 3.784974222639325*^9}, {3.7849745434126034`*^9, 
   3.784974544110735*^9}, {3.784974601321229*^9, 3.784974605461156*^9}, {
   3.784974710736743*^9, 3.7849747684423428`*^9}, {3.784976282298376*^9, 
   3.7849762963519273`*^9}, {3.7849764993822412`*^9, 3.784976515862524*^9}, {
   3.784976673077045*^9, 3.7849766737662954`*^9}, {3.7849770582390857`*^9, 
   3.7849770817050214`*^9}, {3.7849771966792464`*^9, 3.784977200744791*^9}, {
   3.7849774477752275`*^9, 3.784977448371665*^9}, {3.784986312824957*^9, 
   3.7849863423274183`*^9}, {3.784987284149398*^9, 3.7849873444393*^9}, {
   3.7849874406239834`*^9, 3.7849874921625605`*^9}, {3.78498791163412*^9, 
   3.784987918550289*^9}, {3.784987956913732*^9, 3.784987957430295*^9}, {
   3.78498830246691*^9, 3.7849883029525776`*^9}, {3.7849884071614265`*^9, 
   3.7849884119066844`*^9}, {3.784988447079715*^9, 3.7849884680815525`*^9}, {
   3.784988533491603*^9, 3.7849885824401703`*^9}, 3.7849886397804585`*^9, {
   3.7850153211825876`*^9, 3.7850153216874065`*^9}, {3.7850156939006624`*^9, 
   3.785015694648699*^9}, 3.7850157255450363`*^9, {3.785015966192396*^9, 
   3.78501597131769*^9}, {3.7850160275173893`*^9, 3.7850160280202823`*^9}, {
   3.785017368432299*^9, 3.78501739604052*^9}, {3.7850174472106466`*^9, 
   3.785017478921421*^9}, {3.7850179338928843`*^9, 3.7850179695007367`*^9}, {
   3.7850180051071196`*^9, 3.7850180061283536`*^9}, {3.785018037973852*^9, 
   3.7850180386849523`*^9}},
 CellLabel->"In[69]:=",ExpressionUUID->"03cf18a0-f1e3-478c-af1d-32b85ea495ba"],

Cell[BoxData[
 RowBox[{"B", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "2", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"4", ",", "2", ",", "2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2", ",", "4", ",", "1"}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.785015255962994*^9, 3.785015301672872*^9}, {
  3.785015420260889*^9, 3.785015437697337*^9}},
 CellLabel->"In[36]:=",ExpressionUUID->"9a0bdda5-fff6-4233-8d25-85399e077918"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B", "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.785017227935514*^9, 3.785017230985494*^9}},
 CellLabel->"In[58]:=",ExpressionUUID->"ab511fbc-18fa-4bdc-acea-d91f643b0124"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "2", "1"},
     {"4", "2", "2"},
     {"2", "4", "1"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.7850177186839848`*^9},
 CellLabel->
  "Out[58]//MatrixForm=",ExpressionUUID->"b26c7702-0cba-40fb-b94c-\
72b1b8c3a1a7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["B"], "Input",
 CellChangeTimes->{3.7850175926359797`*^9},
 CellLabel->"In[57]:=",ExpressionUUID->"c2fc334d-9952-414a-ba35-3b4dfd2b27f9"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "2", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4", ",", "2", ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2", ",", "4", ",", "1"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.7850152602108474`*^9, 3.7850153025466824`*^9, 
  3.785015438939251*^9, 3.785017231489148*^9, 3.785017594027297*^9},
 CellLabel->"Out[57]=",ExpressionUUID->"7b46d1ba-0090-4ece-8357-9e57b7b2b8a5"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"GuassJordon", "[", "B", "]"}], "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.7850152375668435`*^9, 3.785015243083733*^9}, {
  3.7850153085227084`*^9, 3.7850153106557283`*^9}},
 CellLabel->"In[70]:=",ExpressionUUID->"51d13254-0df7-4704-bf49-d6286f85ad74"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0"},
     {"0", "1", "0"},
     {"0", "0", "1"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.7850152442572856`*^9, {3.785015305838398*^9, 3.785015326846175*^9}, 
   3.7850154433903637`*^9, {3.7850157041945057`*^9, 3.7850157305039477`*^9}, 
   3.78501598095651*^9, 3.7850160321605873`*^9, 3.7850162776235013`*^9, 
   3.7850171712243776`*^9, 3.7850175051235957`*^9, {3.78501799163955*^9, 
   3.7850180145478587`*^9}, 3.785018046190499*^9},
 CellLabel->
  "Out[70]//MatrixForm=",ExpressionUUID->"a0f1a3bf-4e2d-41f0-9bba-\
91e7304a71bd"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B", "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.785017788626028*^9, 3.7850177916299987`*^9}},
 CellLabel->"In[60]:=",ExpressionUUID->"663a59d5-690f-4c7d-8e78-d9e8a1a18fec"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "2", "1"},
     {"4", "2", "2"},
     {"2", "4", "1"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.7850177920548673`*^9},
 CellLabel->
  "Out[60]//MatrixForm=",ExpressionUUID->"897ba509-0957-4b4d-aabe-\
113893b5f538"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"First", "[", 
   RowBox[{
    RowBox[{"Position", "[", 
     RowBox[{
      RowBox[{"Abs", "[", 
       RowBox[{"B", "[", 
        RowBox[{
         RowBox[{"Table", "[", 
          RowBox[{"j", ",", 
           RowBox[{"{", 
            RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", "1"}], 
        "]"}], "]"}], ",", "4"}], "]"}], "//", "Flatten"}], "]"}], "+", 
  RowBox[{"(", 
   RowBox[{"1", "-", "1"}], ")"}]}]], "Input",
 CellChangeTimes->{{3.7850176942082815`*^9, 3.785017728113247*^9}, {
  3.785017848735933*^9, 3.7850178495268526`*^9}},
 CellLabel->"In[61]:=",ExpressionUUID->"734c34da-98a4-4a5a-9511-b66e4e4a216d"],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.785017731750523*^9, 3.7850178507994537`*^9},
 CellLabel->"Out[61]=",ExpressionUUID->"12bb7bd1-88b0-469b-8c2f-679f0613800a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B", "[", 
  RowBox[{
   RowBox[{"Table", "[", 
    RowBox[{"j", ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", "1"}], 
  "]"}]], "Input",
 CellLabel->"In[62]:=",ExpressionUUID->"2b64211b-bf82-4c8f-863d-a458f87781e5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "2", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"4", ",", "2", ",", "2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2", ",", "4", ",", "1"}], "}"}]}], "}"}], "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "2", ",", "3"}], "}"}], ",", "1"}], "]"}]], "Output",
 CellChangeTimes->{3.785017865936103*^9},
 CellLabel->"Out[62]=",ExpressionUUID->"98129c92-b617-497c-8e55-169e764a3ed6"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{"j", ",", 
   RowBox[{"{", 
    RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]], "Input",
 CellLabel->"In[63]:=",ExpressionUUID->"187d9754-897c-4f34-84ae-95937a6d96d5"],

Cell[BoxData[
 RowBox[{"B", "[", 
  RowBox[{"[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "2", ",", "3"}], "}"}], ",", "1"}], "]"}], "]"}]], \
"Input",
 CellChangeTimes->{{3.7850179080724363`*^9, 3.7850179164759803`*^9}},
 CellLabel->"In[64]:=",ExpressionUUID->"2513bbad-9528-41ae-b06c-8ca0a2ec8c7b"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "4", ",", "2"}], "}"}]], "Output",
 CellChangeTimes->{3.7850179201271777`*^9},
 CellLabel->"Out[64]=",ExpressionUUID->"f520807b-dc45-476b-bbbb-79d9c13c130c"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Inverse", "[", "B", "]"}], "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.7850153333633347`*^9, 3.7850153462709694`*^9}},
 CellLabel->"In[35]:=",ExpressionUUID->"8ebc8e24-1a22-44f6-a9ad-b90763b20b01"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", 
       FractionBox["13181", "110875"]}], 
      RowBox[{"-", 
       FractionBox["693", "110875"]}], 
      FractionBox["13682", "110875"], 
      FractionBox["768", "110875"]},
     {
      FractionBox["226", "110875"], 
      RowBox[{"-", 
       FractionBox["47", "110875"]}], 
      RowBox[{"-", 
       FractionBox["672", "110875"]}], 
      FractionBox["1972", "110875"]},
     {
      RowBox[{"-", 
       FractionBox["1583", "110875"]}], 
      FractionBox["1801", "110875"], 
      RowBox[{"-", 
       FractionBox["199", "110875"]}], 
      RowBox[{"-", 
       FractionBox["76", "110875"]}]},
     {
      FractionBox["31636", "110875"], 
      RowBox[{"-", 
       FractionBox["692", "110875"]}], 
      RowBox[{"-", 
       FractionBox["2817", "110875"]}], 
      RowBox[{"-", 
       FractionBox["1633", "110875"]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.7850153385065174`*^9, 3.7850153467044735`*^9}},
 CellLabel->
  "Out[35]//MatrixForm=",ExpressionUUID->"e7a72ae7-a23c-40a2-b4e1-\
47cdee59c48c"]
}, Open  ]],

Cell[" ", "Text",
 Editable->False,
 Selectable->False,
 CellFrame->{{0, 0}, {0, 2}},
 ShowCellBracket->False,
 CellMargins->{{0, 0}, {1, 1}},
 CellElementSpacings->{"CellMinHeight"->1},
 CellFrameMargins->0,
 CellFrameColor->RGBColor[0, 0, 1],
 CellSize->{
  Inherited, 4},ExpressionUUID->"c3562191-a342-44d0-842e-b936de6643f4"],

Cell[CellGroupData[{

Cell["LU\:5206\:89e3\:ff08\:6210\:529f\:ff09", "Section",
 CellChangeTimes->{{3.7849777158171577`*^9, 3.7849777188888245`*^9}, {
  3.7850143212949553`*^9, 
  3.7850143252535205`*^9}},ExpressionUUID->"544df9ba-2d11-4896-b912-\
9d4bb010df07"],

Cell[BoxData[
 RowBox[{
  RowBox[{"LUSolver", "[", 
   RowBox[{"A_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "c", ",", "r", ",", "U", ",", "L", ",", "k", ",", "i", ",", "j", ",", 
      "i1", ",", "i2", ",", "x", ",", "y", ",", "get1", ",", "get2", ",", 
      "n"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"c", "=", 
      RowBox[{"First", "[", 
       RowBox[{"Dimensions", "[", "A", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"U", "=", 
      RowBox[{"UpperTriangularize", "[", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"i", "+", "j"}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "c"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "c"}], "}"}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"L", "=", 
      RowBox[{"IdentityMatrix", "[", "c", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"U", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "All"}], "]"}], "]"}], "=", 
      RowBox[{"A", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "All"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"L", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "1"}], "]"}], "]"}], "=", 
      RowBox[{
       RowBox[{"A", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "1"}], "]"}], "]"}], "/", 
       RowBox[{"U", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "1"}], "]"}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"L", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", "1"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"r", "=", "2"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"r", "<", 
        RowBox[{"c", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"i1", "=", "r"}], ";", 
          RowBox[{"While", "[", 
           RowBox[{
            RowBox[{"i1", "<", "c"}], ",", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{
                RowBox[{"U", "[", 
                 RowBox[{"[", 
                  RowBox[{"r", ",", "i1"}], "]"}], "]"}], "=", 
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"r", ",", "i1"}], "]"}], "]"}], "-", 
                 RowBox[{"Sum", "[", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"L", "[", 
                    RowBox[{"[", 
                    RowBox[{"r", ",", "k"}], "]"}], "]"}], "*", 
                    RowBox[{"U", "[", 
                    RowBox[{"[", 
                    RowBox[{"k", ",", "i1"}], "]"}], "]"}]}], ",", 
                   RowBox[{"{", 
                    RowBox[{"k", ",", "1", ",", 
                    RowBox[{"r", "-", "1"}]}], "}"}]}], "]"}]}]}], ",", 
               RowBox[{
                RowBox[{"i2", "=", 
                 RowBox[{"r", "+", "1"}]}], ";", 
                RowBox[{
                 RowBox[{"L", "[", 
                  RowBox[{"[", 
                   RowBox[{
                    RowBox[{"i1", "+", "1"}], ",", "r"}], "]"}], "]"}], "=", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{
                    RowBox[{"i1", "+", "1"}], ",", "r"}], "]"}], "]"}], "-", 
                    RowBox[{"Sum", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"L", "[", 
                    RowBox[{"[", 
                    RowBox[{
                    RowBox[{"i1", "+", "1"}], ",", "k"}], "]"}], "]"}], "*", 
                    RowBox[{"U", "[", 
                    RowBox[{"[", 
                    RowBox[{"k", ",", "r"}], "]"}], "]"}]}], ",", 
                    RowBox[{"{", 
                    RowBox[{"k", ",", "1", ",", 
                    RowBox[{"r", "-", "1"}]}], "}"}]}], "]"}]}], ")"}], "/", 
                  RowBox[{"U", "[", 
                   RowBox[{"[", 
                    RowBox[{"r", ",", "r"}], "]"}], "]"}]}]}]}]}], "}"}], ";", 
             RowBox[{"i1", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"U", "[", 
            RowBox[{"[", 
             RowBox[{"r", ",", "c"}], "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"A", "[", 
             RowBox[{"[", 
              RowBox[{"r", ",", "c"}], "]"}], "]"}], "-", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"L", "[", 
                RowBox[{"[", 
                 RowBox[{"r", ",", "k"}], "]"}], "]"}], "*", 
               RowBox[{"U", "[", 
                RowBox[{"[", 
                 RowBox[{"k", ",", "c"}], "]"}], "]"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "1", ",", 
                RowBox[{"r", "-", "1"}]}], "}"}]}], "]"}]}]}]}], "}"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"r", "++"}]}]}], "\[IndentingNewLine]", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"y", "=", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"x", "=", "y"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"y", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "=", 
      RowBox[{"b", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"i", "=", "2"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"i", "<", 
        RowBox[{"c", "+", "1"}]}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{"y", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "=", 
         RowBox[{
          RowBox[{"b", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "-", 
          RowBox[{"Sum", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"L", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "k"}], "]"}], "]"}], "*", 
             RowBox[{"y", "[", 
              RowBox[{"[", "k", "]"}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"k", ",", "1", ",", 
              RowBox[{"i", "-", "1"}]}], "}"}]}], "]"}]}]}], ";", 
        RowBox[{"i", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"x", "[", 
       RowBox[{"[", "c", "]"}], "]"}], "=", 
      RowBox[{
       RowBox[{"y", "[", 
        RowBox[{"[", "c", "]"}], "]"}], "/", 
       RowBox[{"U", "[", 
        RowBox[{"[", 
         RowBox[{"c", ",", "c"}], "]"}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"j", "=", 
      RowBox[{"c", "-", "1"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"j", ">", "0"}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{"x", "[", 
          RowBox[{"[", "j", "]"}], "]"}], "=", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"y", "[", 
             RowBox[{"[", "j", "]"}], "]"}], "-", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"U", "[", 
                RowBox[{"[", 
                 RowBox[{"j", ",", "k"}], "]"}], "]"}], "*", 
               RowBox[{"x", "[", 
                RowBox[{"[", "k", "]"}], "]"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", 
                RowBox[{"j", "+", "1"}], ",", "c"}], "}"}]}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"U", "[", 
           RowBox[{"[", 
            RowBox[{"j", ",", "j"}], "]"}], "]"}]}]}], ";", 
        RowBox[{"j", "--"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Return", "[", "x", "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.7849789276688805`*^9, 3.7849789283252487`*^9}, 
   3.784978980908723*^9, {3.784988950395756*^9, 3.784988950985146*^9}, 
   3.785014312262432*^9, {3.7850314639386663`*^9, 3.7850314670532093`*^9}, {
   3.7850316808032746`*^9, 3.7850318885573387`*^9}, {3.7850319354899197`*^9, 
   3.785031988142519*^9}, {3.7850321266485095`*^9, 3.785032173903014*^9}, {
   3.7850322261697283`*^9, 3.785032321798702*^9}, {3.7850328638226395`*^9, 
   3.785032875930019*^9}, {3.7850332018124576`*^9, 3.785033295358678*^9}, {
   3.785033382496922*^9, 3.7850333829591603`*^9}, {3.785033432283819*^9, 
   3.7850334510011344`*^9}, {3.785033586819111*^9, 3.7850338995755925`*^9}, {
   3.785033935044347*^9, 3.785033962623983*^9}, 3.785034014932987*^9, {
   3.785034502692917*^9, 3.7850345031312094`*^9}},
 CellLabel->"In[88]:=",ExpressionUUID->"3d281545-cce5-40a1-9b8c-9d3906fa76e4"],

Cell[CellGroupData[{

Cell[BoxData["B"], "Input",
 CellChangeTimes->{3.7849789008191566`*^9},
 CellLabel->"In[85]:=",ExpressionUUID->"eb3ac092-d10c-40cf-ae97-ef6901f6f0ff"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "6", ",", "21", ",", "666"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"7", ",", "4", ",", "2", ",", "56"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2", ",", "6", ",", "11", ",", "17"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3", ",", "99", ",", "102", ",", "1"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.784978901560132*^9},
 CellLabel->"Out[85]=",ExpressionUUID->"9d475d11-c9cd-41dc-a54d-8ab9b738409a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"LUSolver", "[", 
   RowBox[{"B", ",", 
    RowBox[{"{", 
     RowBox[{"1", ",", "2", ",", "3"}], "}"}]}], "]"}], "//", 
  "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.784978907207222*^9, 3.784978914771533*^9}, {
  3.7850333124280734`*^9, 3.785033315786721*^9}, {3.7850334611901793`*^9, 
  3.7850334639807577`*^9}},
 CellLabel->"In[89]:=",ExpressionUUID->"5d92213f-562b-4370-b4b9-c86987497bd4"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {
       FractionBox["2", "3"]},
      {
       FractionBox["2", "3"]},
      {
       RowBox[{"-", "1"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.784978915198186*^9, 3.7849789328408384`*^9}, 
   3.78497902385868*^9, 3.7850333179458766`*^9, 3.7850333881449184`*^9, {
   3.7850334371287985`*^9, 3.7850334663224335`*^9}, 3.7850340488000684`*^9, 
   3.7850345100844946`*^9},
 CellLabel->
  "Out[89]//MatrixForm=",ExpressionUUID->"116b9dae-6cad-4a0d-8ce0-\
62b664e6f91f"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["B"], "Input",
 CellChangeTimes->{3.78503339486534*^9},
 CellLabel->"In[79]:=",ExpressionUUID->"a953f32c-a9d1-4965-b2e4-be2bea9a1512"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "2", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4", ",", "2", ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2", ",", "4", ",", "1"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.785033396872815*^9},
 CellLabel->"Out[79]=",ExpressionUUID->"0f144c48-9519-4e5c-9345-213c7cf6ef57"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"LUSolve", "[", "A_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "c", ",", "r", ",", "U", ",", "L", ",", "k", ",", "i", ",", "i1", ",", 
      "i2", ",", "x", ",", "y"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"c", "=", 
      RowBox[{"First", "[", 
       RowBox[{"Dimensions", "[", "A", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"U", "=", 
      RowBox[{"UpperTriangularize", "[", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"i", "+", "j"}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "c"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "c"}], "}"}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"L", "=", 
      RowBox[{"IdentityMatrix", "[", "c", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"U", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "All"}], "]"}], "]"}], "=", 
      RowBox[{"A", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "All"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"L", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "1"}], "]"}], "]"}], "=", 
      RowBox[{
       RowBox[{"A", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "1"}], "]"}], "]"}], "/", 
       RowBox[{"U", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "1"}], "]"}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"L", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", "1"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"r", "=", "2"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"r", "<", 
        RowBox[{"c", "+", "1"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"i1", "=", "r"}], ";", 
          RowBox[{"While", "[", 
           RowBox[{
            RowBox[{"i1", "<", 
             RowBox[{"c", "+", "1"}]}], ",", 
            RowBox[{
             RowBox[{
              RowBox[{"U", "[", 
               RowBox[{"[", 
                RowBox[{"r", ",", "i1"}], "]"}], "]"}], "=", 
              RowBox[{
               RowBox[{"A", "[", 
                RowBox[{"[", 
                 RowBox[{"r", ",", "i1"}], "]"}], "]"}], "-", 
               RowBox[{"Sum", "[", 
                RowBox[{
                 RowBox[{
                  RowBox[{"L", "[", 
                   RowBox[{"[", 
                    RowBox[{"r", ",", "k"}], "]"}], "]"}], "*", 
                  RowBox[{"U", "[", 
                   RowBox[{"[", 
                    RowBox[{"k", ",", "i1"}], "]"}], "]"}]}], ",", 
                 RowBox[{"{", 
                  RowBox[{"k", ",", "1", ",", 
                   RowBox[{"r", "-", "1"}]}], "}"}]}], "]"}]}]}], ";", 
             RowBox[{"i1", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
          RowBox[{"i2", "=", 
           RowBox[{"r", "+", "1"}]}], ";", 
          RowBox[{"While", "[", 
           RowBox[{
            RowBox[{"i2", "<", 
             RowBox[{"c", "+", "1"}]}], ",", 
            RowBox[{
             RowBox[{
              RowBox[{"L", "[", 
               RowBox[{"[", 
                RowBox[{"i2", ",", "r"}], "]"}], "]"}], "=", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"A", "[", 
                  RowBox[{"[", 
                   RowBox[{"i2", ",", "r"}], "]"}], "]"}], "-", 
                 RowBox[{"Sum", "[", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"L", "[", 
                    RowBox[{"[", 
                    RowBox[{"i2", ",", "k"}], "]"}], "]"}], "*", 
                    RowBox[{"U", "[", 
                    RowBox[{"[", 
                    RowBox[{"k", ",", "r"}], "]"}], "]"}]}], ",", 
                   RowBox[{"{", 
                    RowBox[{"k", ",", "1", ",", 
                    RowBox[{"r", "-", "1"}]}], "}"}]}], "]"}]}], ")"}], "/", 
               RowBox[{"U", "[", 
                RowBox[{"[", 
                 RowBox[{"r", ",", "r"}], "]"}], "]"}]}]}], ";", 
             RowBox[{"i2", "++"}]}]}], "]"}]}], "}"}], ";", 
        RowBox[{"r", "++"}]}]}], "\[IndentingNewLine]", "\[IndentingNewLine]",
       "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"y", "=", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"x", "=", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"y", "[", 
      RowBox[{"[", "i", "]"}], "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.7850314155294323`*^9, 3.785031439321149*^9}, {
  3.785031485027963*^9, 3.7850314886595793`*^9}, {3.7850315243223867`*^9, 
  3.7850316369756913`*^9}},ExpressionUUID->"44bdc003-a7c7-4e48-ba61-\
958730091d5a"],

Cell[" ", "Text",
 Editable->False,
 Selectable->False,
 CellFrame->{{0, 0}, {0, 2}},
 ShowCellBracket->False,
 CellMargins->{{0, 0}, {1, 1}},
 CellElementSpacings->{"CellMinHeight"->1},
 CellFrameMargins->0,
 CellFrameColor->RGBColor[0, 0, 1],
 CellSize->{
  Inherited, 4},ExpressionUUID->"931b6f6b-47e3-4ed2-b7f2-f6c41a4667e7"],

Cell[" ", "Text",
 Editable->False,
 Selectable->False,
 CellFrame->{{0, 0}, {0, 2}},
 ShowCellBracket->False,
 CellMargins->{{0, 0}, {1, 1}},
 CellElementSpacings->{"CellMinHeight"->1},
 CellFrameMargins->0,
 CellFrameColor->RGBColor[0, 0, 1],
 CellSize->{
  Inherited, 4},ExpressionUUID->"6db80a03-209f-42e1-9d49-8435f4a92963"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\:8fed\:4ee3\:6cd5", "Section",
 CellChangeTimes->{{3.7849798719556875`*^9, 
  3.784979873504547*^9}},ExpressionUUID->"ab435dfa-bde9-4bf5-975d-\
3e57cf37e8da"],

Cell[CellGroupData[{

Cell["\:96c5\:514b\:6bd4\:8fed\:4ee3\:6cd5\:ff08\:6210\:529f\:ff09", \
"Subsection",
 CellChangeTimes->{{3.7849798815470486`*^9, 3.7849798840244284`*^9}, {
  3.784984699065511*^9, 
  3.7849847044531093`*^9}},ExpressionUUID->"0398b880-538e-4f2e-8967-\
143e7eec2da6"],

Cell[BoxData[
 RowBox[{
  RowBox[{"JacobiNest", "[", 
   RowBox[{"A_", ",", "b_", ",", "rand_", ",", "ops_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"B", "=", 
       RowBox[{
        RowBox[{"-", "A"}], "+", 
        RowBox[{"IdentityMatrix", "[", 
         RowBox[{"Dimensions", "[", "A", "]"}], "]"}]}]}], ",", 
      RowBox[{"f", "=", "b"}], ",", "y"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"y", "=", 
      RowBox[{"NestWhile", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"B", ".", "#"}], "+", "f"}], ")"}], "&"}], ",", "rand", 
        ",", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Abs", "[", 
            RowBox[{
             RowBox[{"Norm", "[", "#1", "]"}], "-", 
             RowBox[{"Norm", "[", "#2", "]"}]}], "]"}], "<", "ops"}], ")"}], 
         "&"}], ",", "2"}], "]"}]}], ";", 
     RowBox[{"Return", "[", "y", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7849799134268684`*^9, 3.78497997247983*^9}, {
  3.7849800347872915`*^9, 3.7849803472235804`*^9}, {3.784980388215968*^9, 
  3.784980395033868*^9}, {3.7849845371728234`*^9, 3.784984540623559*^9}},
 CellLabel->
  "In[134]:=",ExpressionUUID->"01c3cd34-de03-42fa-a4a2-72ef6acc146b"],

Cell[" ", "Text",
 Editable->False,
 Selectable->False,
 CellFrame->{{0, 0}, {0, 2}},
 ShowCellBracket->False,
 CellMargins->{{0, 0}, {1, 1}},
 CellElementSpacings->{"CellMinHeight"->1},
 CellFrameMargins->0,
 CellFrameColor->RGBColor[0, 0, 1],
 CellSize->{
  Inherited, 4},ExpressionUUID->"4918d143-afdc-472a-924b-c798cf3f2ccd"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\:9ad8\:65af-\:8d5b\:5fb7\:5c14\:8fed\:4ee3\:6cd5\:ff08\:6210\:529f\
\:ff09", "Subsection",
 CellChangeTimes->{{3.784991822967889*^9, 3.784991831836345*^9}, {
  3.7850099408116064`*^9, 
  3.7850099457892995`*^9}},ExpressionUUID->"10b3b4bc-2a60-419f-8c69-\
026bc319969d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"GuassSeidle", "[", 
   RowBox[{"A_", ",", "b_", ",", "rand_", ",", "ops_"}], "]"}], ":=", 
  "\[IndentingNewLine]", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"c", "=", 
       RowBox[{"First", "[", 
        RowBox[{"Dimensions", "[", "A", "]"}], "]"}]}], ",", "x2", ",", "x1", 
      ",", "n"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"x1", "=", "rand"}], ";", "\[IndentingNewLine]", 
     RowBox[{"n", "=", "1"}], ";", "\[IndentingNewLine]", 
     RowBox[{"x2", "=", "rand"}], ";", "\[IndentingNewLine]", 
     "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Abs", "[", 
         RowBox[{
          RowBox[{"Norm", "[", "x1", "]"}], "-", 
          RowBox[{"Norm", "[", "x2", "]"}]}], "]"}], "<", "ops"}], ",", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"x2", "=", 
           RowBox[{"Table", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"b", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "-", 
                RowBox[{"Sum", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
                   RowBox[{"x2", "[", 
                    RowBox[{"[", "j", "]"}], "]"}]}], ",", 
                  RowBox[{"{", 
                   RowBox[{"j", ",", "1", ",", 
                    RowBox[{"i", "-", "1"}]}], "}"}]}], "]"}], "-", 
                RowBox[{"Sum", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"A", "[", 
                    RowBox[{"[", 
                    RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
                   RowBox[{"x1", "[", 
                    RowBox[{"[", "j", "]"}], "]"}]}], ",", 
                  RowBox[{"{", 
                   RowBox[{"j", ",", 
                    RowBox[{"i", "+", "1"}], ",", "c"}], "}"}]}], "]"}]}], 
               ")"}], "/", 
              RowBox[{"A", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "i"}], "]"}], "]"}]}], ",", 
             RowBox[{"{", 
              RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"x1", ",", "x2"}], "}"}], "=", 
           RowBox[{"{", 
            RowBox[{"x2", ",", "x1"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
         "}"}], ";", "\[IndentingNewLine]", 
        RowBox[{"n", "++"}]}]}], "\[IndentingNewLine]", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Return", "[", "x1", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.784991845763427*^9, 3.78499210056188*^9}, 
   3.7849936327785997`*^9, {3.7849936767101264`*^9, 3.7849940433122053`*^9}, {
   3.7849940789299912`*^9, 3.784994245581234*^9}, {3.784994304625401*^9, 
   3.784994319449773*^9}, {3.7849943608036065`*^9, 3.7849943610339923`*^9}, {
   3.784994457217993*^9, 3.784994561675762*^9}, {3.784994591758393*^9, 
   3.7849945988992605`*^9}, {3.7849947031750603`*^9, 
   3.7849947060134335`*^9}, {3.785008835704487*^9, 3.7850089038387175`*^9}, {
   3.7850089817321806`*^9, 3.7850089832640824`*^9}, 3.7850090403276587`*^9, {
   3.7850093197880754`*^9, 3.785009339985087*^9}, {3.785009924773965*^9, 
   3.785009925585795*^9}},
 CellLabel->"In[24]:=",ExpressionUUID->"99913f9a-9c8f-4532-8019-26c7c1448f3a"],

Cell[" ", "Text",
 Editable->False,
 Selectable->False,
 CellFrame->{{0, 0}, {0, 2}},
 ShowCellBracket->False,
 CellMargins->{{0, 0}, {1, 1}},
 CellElementSpacings->{"CellMinHeight"->1},
 CellFrameMargins->0,
 CellFrameColor->RGBColor[0, 0, 1],
 CellSize->{
  Inherited, 4},ExpressionUUID->"3ae9cd2e-261e-4529-9615-c6ae4663ea64"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\:8ffd\:8d76\:6cd5", "Section",
 CellChangeTimes->{{3.785035219847252*^9, 
  3.7850352269011607`*^9}},ExpressionUUID->"1c2529f6-eba0-49be-a041-\
cf15b9ec240f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Chasing", "[", 
   RowBox[{"A_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "c", ",", "d", ",", "u", ",", "l", ",", "i", ",", "n", ",", "x", ",", 
      "k"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"c", "=", 
      RowBox[{
       RowBox[{
        RowBox[{"Dimensions", "[", "A", "]"}], "//", "Flatten"}], "//", 
       "First"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"d", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"A", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "i"}], "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"u", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"A", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "-", "1"}], ",", "i"}], "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "2", ",", "c"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"l", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"A", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", 
           RowBox[{"i", "+", "1"}]}], "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", 
          RowBox[{"c", "-", "1"}]}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"i", "=", "2"}], ";", "\[IndentingNewLine]", 
     RowBox[{"k", "=", "b"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"i", "<", 
        RowBox[{"c", "+", "1"}]}], ",", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"l", "[", 
            RowBox[{"[", 
             RowBox[{"i", "-", "1"}], "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"l", "[", 
             RowBox[{"[", 
              RowBox[{"i", "-", "1"}], "]"}], "]"}], "/", 
            RowBox[{"d", "[", 
             RowBox[{"[", 
              RowBox[{"i", "-", "1"}], "]"}], "]"}]}]}], ";", 
          RowBox[{
           RowBox[{"d", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"d", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "-", 
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", 
               RowBox[{"i", "-", "1"}], "]"}], "]"}], "*", 
             RowBox[{"l", "[", 
              RowBox[{"[", 
               RowBox[{"i", "-", "1"}], "]"}], "]"}]}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"k", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "=", 
           RowBox[{
            RowBox[{"k", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "-", 
            RowBox[{
             RowBox[{"k", "[", 
              RowBox[{"[", 
               RowBox[{"i", "-", "1"}], "]"}], "]"}], "*", 
             RowBox[{"l", "[", 
              RowBox[{"[", 
               RowBox[{"i", "-", "1"}], "]"}], "]"}]}]}]}]}], "}"}], ";", 
        RowBox[{"i", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"n", "=", 
      RowBox[{"c", "-", "1"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"x", "=", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"x", "[", 
       RowBox[{"[", "c", "]"}], "]"}], "=", 
      RowBox[{
       RowBox[{"k", "[", 
        RowBox[{"[", "c", "]"}], "]"}], "/", 
       RowBox[{"d", "[", 
        RowBox[{"[", "c", "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{"n", ">", "0"}], ",", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"x", "[", 
           RowBox[{"[", "n", "]"}], "]"}], "=", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"k", "[", 
              RowBox[{"[", "n", "]"}], "]"}], "-", 
             RowBox[{
              RowBox[{"u", "[", 
               RowBox[{"[", "n", "]"}], "]"}], "*", 
              RowBox[{"x", "[", 
               RowBox[{"[", 
                RowBox[{"n", "+", "1"}], "]"}], "]"}]}]}], ")"}], "/", 
           RowBox[{"d", "[", 
            RowBox[{"[", "n", "]"}], "]"}]}]}], "}"}], ";", 
        RowBox[{"n", "--"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Return", "[", "x", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7850352338204427`*^9, 3.7850352512777767`*^9}, {
  3.785035287794695*^9, 3.7850358024257593`*^9}, {3.785036848192683*^9, 
  3.785036868976287*^9}, {3.785037696716958*^9, 3.785037699475154*^9}, {
  3.785038086084979*^9, 3.7850381031711264`*^9}, {3.7850391196885867`*^9, 
  3.785039134283274*^9}, {3.7850394459895096`*^9, 3.785039448958105*^9}, {
  3.785039498795226*^9, 3.785039503977104*^9}, {3.7850395474506006`*^9, 
  3.7850395492984295`*^9}, {3.7850397936671953`*^9, 3.785039798911227*^9}, {
  3.785039833314088*^9, 3.785039836810652*^9}, {3.7850398670527477`*^9, 
  3.7850398674228306`*^9}, {3.7850399260041623`*^9, 3.7850399633646154`*^9}},
 CellLabel->
  "In[169]:=",ExpressionUUID->"9c2dd963-5d83-495a-8b25-173af639ecfa"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"s", "=", 
  RowBox[{"SparseArray", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"i_", ",", "i_"}], "}"}], "\[Rule]", 
       RowBox[{"-", "2"}]}], ",", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{"i_", ",", "j_"}], "}"}], "/;", 
        RowBox[{
         RowBox[{"Abs", "[", 
          RowBox[{"i", "-", "j"}], "]"}], "\[Equal]", "1"}]}], "\[Rule]", 
       "1"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"5", ",", "5"}], "}"}]}], "]"}]}]], "Input",
 CellLabel->"In[91]:=",ExpressionUUID->"03fcadc4-b839-4aeb-be71-dd81d9e7a5b5"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   TagBox["SparseArray",
    "SummaryHead"], "[", 
   DynamicModuleBox[{Typeset`open$$ = False, Typeset`embedState$$ = "Ready"}, 
    TemplateBox[{PaneSelectorBox[{False -> GridBox[{{
            PaneBox[
             ButtonBox[
              DynamicBox[
               FEPrivate`FrontEndResource[
               "FEBitmaps", "SquarePlusIconMedium"]], 
              ButtonFunction :> (Typeset`open$$ = True), Appearance -> None, 
              Evaluator -> Automatic, Method -> "Preemptive"], 
             Alignment -> {Center, Center}, ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}]], 
            GraphicsBox[
             RasterBox[CompressedData["
1:eJxTTMoPSmVmYGBghWIQ26p8d9wspXP2c4wXbTH/ccX+xSLFx4GrXtmTKj6/
yEMzreG+vaNe4ZPGEy/sz/2orS2ofW+v8D5UW7juqD0DFIhsOflwst9lmpuP
S/1wNR+XOK3NxyUOAC4TJC4=
              "], {{0, 0}, {5, 5}}, {0, 1}], AspectRatio -> 1, Axes -> False, 
             Frame -> False, FrameLabel -> {None, None}, FrameStyle -> 
             Directive[
               Opacity[0.5], 
               Thickness[Tiny], 
               RGBColor[0.368417, 0.506779, 0.709798]], 
             FrameTicks -> {{None, None}, {None, None}}, GridLinesStyle -> 
             Directive[
               GrayLevel[0.5, 0.4]], ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}], 
             Method -> {
              "AxisPadding" -> Scaled[0.02], "DefaultBoundaryStyle" -> 
               Automatic, "DefaultPlotStyle" -> Automatic, "DomainPadding" -> 
               Scaled[0.02], "RangePadding" -> Scaled[0.05]}], 
            GridBox[{{
               RowBox[{
                 TagBox["\"Specified elements: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["13", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Dimensions: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox[
                  RowBox[{"{", 
                    RowBox[{"5", ",", "5"}], "}"}], "SummaryItem"]}]}}, 
             GridBoxAlignment -> {
              "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> 
             False, GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
              BaseStyle -> {
              ShowStringCharacters -> False, NumberMarks -> False, 
               PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, 
          GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, 
          GridBoxItemSize -> {
           "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
          BaselinePosition -> {1, 1}], True -> GridBox[{{
            PaneBox[
             ButtonBox[
              DynamicBox[
               FEPrivate`FrontEndResource[
               "FEBitmaps", "SquareMinusIconMedium"]], 
              ButtonFunction :> (Typeset`open$$ = False), Appearance -> None, 
              Evaluator -> Automatic, Method -> "Preemptive"], 
             Alignment -> {Center, Center}, ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}]], 
            GraphicsBox[
             RasterBox[CompressedData["
1:eJxTTMoPSmVmYGBghWIQ26p8d9wspXP2c4wXbTH/ccX+xSLFx4GrXtmTKj6/
yEMzreG+vaNe4ZPGEy/sz/2orS2ofW+v8D5UW7juqD0DFIhsOflwst9lmpuP
S/1wNR+XOK3NxyUOAC4TJC4=
              "], {{0, 0}, {5, 5}}, {0, 1}], AspectRatio -> 1, Axes -> False, 
             Frame -> False, FrameLabel -> {None, None}, FrameStyle -> 
             Directive[
               Opacity[0.5], 
               Thickness[Tiny], 
               RGBColor[0.368417, 0.506779, 0.709798]], 
             FrameTicks -> {{None, None}, {None, None}}, GridLinesStyle -> 
             Directive[
               GrayLevel[0.5, 0.4]], ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}], 
             Method -> {
              "AxisPadding" -> Scaled[0.02], "DefaultBoundaryStyle" -> 
               Automatic, "DefaultPlotStyle" -> Automatic, "DomainPadding" -> 
               Scaled[0.02], "RangePadding" -> Scaled[0.05]}], 
            GridBox[{{
               RowBox[{
                 TagBox["\"Specified elements: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["13", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Dimensions: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox[
                  RowBox[{"{", 
                    RowBox[{"5", ",", "5"}], "}"}], "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Default: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["0", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Density: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["0.52`", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Elements:\"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["\"\"", "SummaryItem"]}]}, {
               TagBox[
                TagBox[
                 GridBox[{{
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"1", ",", "2"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    "1", ContentPadding -> False, FrameMargins -> 0, 
                    StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"1", ",", "1"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    RowBox[{"-", "2"}], ContentPadding -> False, FrameMargins -> 
                    0, StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"2", ",", "1"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    "1", ContentPadding -> False, FrameMargins -> 0, 
                    StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"2", ",", "2"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    RowBox[{"-", "2"}], ContentPadding -> False, FrameMargins -> 
                    0, StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    ItemBox[
                    "\"\[VerticalEllipsis]\"", Alignment -> Center, 
                    StripOnInput -> False]}}, 
                  GridBoxAlignment -> {"Columns" -> {{Left}}}, 
                  DefaultBaseStyle -> "Column", 
                  GridBoxItemSize -> {
                   "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], 
                 "Column"], "SummaryItem"]}}, 
             GridBoxAlignment -> {
              "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> 
             False, GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
              BaseStyle -> {
              ShowStringCharacters -> False, NumberMarks -> False, 
               PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, 
          GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, 
          GridBoxItemSize -> {
           "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
          BaselinePosition -> {1, 1}]}, 
       Dynamic[Typeset`open$$], ImageSize -> Automatic]},
     "SummaryPanel"],
    DynamicModuleValues:>{}], "]"}],
  SparseArray[
  Automatic, {5, 5}, 0, {
   1, {{0, 2, 5, 8, 11, 13}, {{2}, {1}, {1}, {2}, {3}, {4}, {3}, {2}, {3}, {
     5}, {4}, {4}, {5}}}, {1, -2, 1, -2, 1, 1, -2, 1, 1, 1, -2, 1, -2}}],
  Editable->False,
  SelectWithContents->True,
  Selectable->False]], "Output",
 CellChangeTimes->{3.785036814394027*^9},
 CellLabel->"Out[91]=",ExpressionUUID->"1cbcbc5e-96c0-4c32-a8e5-63a0e47971ad"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Dimensions", "[", "s", "]"}], "//", "Flatten"}], "//", 
  "First"}]], "Input",
 CellChangeTimes->{{3.7850369460769653`*^9, 3.785036954829285*^9}},
 CellLabel->"In[95]:=",ExpressionUUID->"3d552d3b-61f1-487d-b536-b3c8bf549289"],

Cell[CellGroupData[{

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.7850394101595435`*^9, 
  3.7850394123531837`*^9}},ExpressionUUID->"20d58777-a045-42db-9dcb-\
e7d954b5f39b"],

Cell[BoxData["5"], "Output",
 CellChangeTimes->{3.7850369554504614`*^9},
 CellLabel->"Out[95]=",ExpressionUUID->"b6d4b11d-a8ea-45a1-abd0-aa4be3eaeb95"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{"1", ",", "2"}], "}"}], "[", 
  RowBox[{"[", "2", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.7850396800921326`*^9, 3.785039692088197*^9}},
 CellLabel->
  "In[163]:=",ExpressionUUID->"81b5ab71-30a3-4ea0-9925-12ca72280e40"],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{{3.7850396872092896`*^9, 3.785039692765628*^9}},
 CellLabel->
  "Out[163]=",ExpressionUUID->"abb9fca9-028a-4a19-8e0b-8912d0fd8521"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LUSolver", "[", 
  RowBox[{"s", ",", 
   RowBox[{"{", 
    RowBox[{"1", ",", "2", ",", "34", ",", "5", ",", "6"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.785040010262336*^9, 3.7850400423984575`*^9}},
 CellLabel->
  "In[173]:=",ExpressionUUID->"b2987884-cf18-49fd-8c43-a84c7dc5f765"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox["131", "6"]}], ",", 
   RowBox[{"-", 
    FractionBox["128", "3"]}], ",", 
   RowBox[{"-", 
    FractionBox["123", "2"]}], ",", 
   RowBox[{"-", 
    FractionBox["139", "3"]}], ",", 
   RowBox[{"-", 
    FractionBox["157", "6"]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.7850400172318373`*^9, 3.7850400437950945`*^9}},
 CellLabel->
  "Out[173]=",ExpressionUUID->"61608589-1733-4257-b2e0-468b406b31c7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Chasing", "[", 
  RowBox[{"s", ",", 
   RowBox[{"{", 
    RowBox[{"1", ",", "2", ",", "34", ",", "5", ",", "6"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.785036818098814*^9, 3.7850368295471077`*^9}},
 CellLabel->
  "In[171]:=",ExpressionUUID->"5a6c3f87-f334-4d9f-be45-aa6e1747f48b"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox["131", "6"]}], ",", 
   RowBox[{"-", 
    FractionBox["128", "3"]}], ",", 
   RowBox[{"-", 
    FractionBox["123", "2"]}], ",", 
   RowBox[{"-", 
    FractionBox["139", "3"]}], ",", 
   RowBox[{"-", 
    FractionBox["157", "6"]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.785036833284067*^9, 3.78503687993412*^9, {3.785037321591215*^9, 
   3.7850373505223475`*^9}, {3.7850377244891987`*^9, 3.7850377310273237`*^9}, 
   3.78503795641147*^9, 3.785038110159867*^9, {3.785039375997774*^9, 
   3.785039385218936*^9}, 3.7850394633566165`*^9, 3.7850395110598545`*^9, 
   3.7850395602486987`*^9, 3.7850398092536907`*^9, 3.7850398437467203`*^9, {
   3.785039980154029*^9, 3.7850399854839563`*^9}},
 CellLabel->
  "Out[171]=",ExpressionUUID->"3775de8f-cad6-4103-8547-2041f887d7bb"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"A", "=", "s"}]], "Input",
 CellChangeTimes->{{3.785037835778237*^9, 3.7850378367676115`*^9}},
 CellLabel->
  "In[103]:=",ExpressionUUID->"5df2f2be-d24c-4e19-ab6b-b64e1d9914e9"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   TagBox["SparseArray",
    "SummaryHead"], "[", 
   DynamicModuleBox[{Typeset`open$$ = False, Typeset`embedState$$ = "Ready"}, 
    TemplateBox[{PaneSelectorBox[{False -> GridBox[{{
            PaneBox[
             ButtonBox[
              DynamicBox[
               FEPrivate`FrontEndResource[
               "FEBitmaps", "SquarePlusIconMedium"]], 
              ButtonFunction :> (Typeset`open$$ = True), Appearance -> None, 
              Evaluator -> Automatic, Method -> "Preemptive"], 
             Alignment -> {Center, Center}, ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}]], 
            GraphicsBox[
             RasterBox[CompressedData["
1:eJxTTMoPSmVmYGBghWIQ26p8d9wspXP2c4wXbTH/ccX+xSLFx4GrXtmTKj6/
yEMzreG+vaNe4ZPGEy/sz/2orS2ofW+v8D5UW7juqD0DFIhsOflwst9lmpuP
S/1wNR+XOK3NxyUOAC4TJC4=
              "], {{0, 0}, {5, 5}}, {0, 1}], AspectRatio -> 1, Axes -> False, 
             Frame -> False, FrameLabel -> {None, None}, FrameStyle -> 
             Directive[
               Opacity[0.5], 
               Thickness[Tiny], 
               RGBColor[0.368417, 0.506779, 0.709798]], 
             FrameTicks -> {{None, None}, {None, None}}, GridLinesStyle -> 
             Directive[
               GrayLevel[0.5, 0.4]], ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}], 
             Method -> {
              "AxisPadding" -> Scaled[0.02], "DefaultBoundaryStyle" -> 
               Automatic, "DefaultPlotStyle" -> Automatic, "DomainPadding" -> 
               Scaled[0.02], "RangePadding" -> Scaled[0.05]}], 
            GridBox[{{
               RowBox[{
                 TagBox["\"Specified elements: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["13", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Dimensions: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox[
                  RowBox[{"{", 
                    RowBox[{"5", ",", "5"}], "}"}], "SummaryItem"]}]}}, 
             GridBoxAlignment -> {
              "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> 
             False, GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
              BaseStyle -> {
              ShowStringCharacters -> False, NumberMarks -> False, 
               PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, 
          GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, 
          GridBoxItemSize -> {
           "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
          BaselinePosition -> {1, 1}], True -> GridBox[{{
            PaneBox[
             ButtonBox[
              DynamicBox[
               FEPrivate`FrontEndResource[
               "FEBitmaps", "SquareMinusIconMedium"]], 
              ButtonFunction :> (Typeset`open$$ = False), Appearance -> None, 
              Evaluator -> Automatic, Method -> "Preemptive"], 
             Alignment -> {Center, Center}, ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}]], 
            GraphicsBox[
             RasterBox[CompressedData["
1:eJxTTMoPSmVmYGBghWIQ26p8d9wspXP2c4wXbTH/ccX+xSLFx4GrXtmTKj6/
yEMzreG+vaNe4ZPGEy/sz/2orS2ofW+v8D5UW7juqD0DFIhsOflwst9lmpuP
S/1wNR+XOK3NxyUOAC4TJC4=
              "], {{0, 0}, {5, 5}}, {0, 1}], AspectRatio -> 1, Axes -> False, 
             Frame -> False, FrameLabel -> {None, None}, FrameStyle -> 
             Directive[
               Opacity[0.5], 
               Thickness[Tiny], 
               RGBColor[0.368417, 0.506779, 0.709798]], 
             FrameTicks -> {{None, None}, {None, None}}, GridLinesStyle -> 
             Directive[
               GrayLevel[0.5, 0.4]], ImageSize -> 
             Dynamic[{
               Automatic, 
                3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                 Magnification])}], 
             Method -> {
              "AxisPadding" -> Scaled[0.02], "DefaultBoundaryStyle" -> 
               Automatic, "DefaultPlotStyle" -> Automatic, "DomainPadding" -> 
               Scaled[0.02], "RangePadding" -> Scaled[0.05]}], 
            GridBox[{{
               RowBox[{
                 TagBox["\"Specified elements: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["13", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Dimensions: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox[
                  RowBox[{"{", 
                    RowBox[{"5", ",", "5"}], "}"}], "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Default: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["0", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Density: \"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["0.52`", "SummaryItem"]}]}, {
               RowBox[{
                 TagBox["\"Elements:\"", "SummaryItemAnnotation"], 
                 "\[InvisibleSpace]", 
                 TagBox["\"\"", "SummaryItem"]}]}, {
               TagBox[
                TagBox[
                 GridBox[{{
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"1", ",", "2"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    "1", ContentPadding -> False, FrameMargins -> 0, 
                    StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"1", ",", "1"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    RowBox[{"-", "2"}], ContentPadding -> False, FrameMargins -> 
                    0, StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"2", ",", "1"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    "1", ContentPadding -> False, FrameMargins -> 0, 
                    StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"2", ",", "2"}], "}"}], "\[Rule]", 
                    StyleBox[
                    PaneBox[
                    RowBox[{"-", "2"}], ContentPadding -> False, FrameMargins -> 
                    0, StripOnInput -> True, BaselinePosition -> Baseline, 
                    ImageSize -> {{1, 300}, Automatic}], LineBreakWithin -> 
                    False]}]}, {
                    ItemBox[
                    "\"\[VerticalEllipsis]\"", Alignment -> Center, 
                    StripOnInput -> False]}}, 
                  GridBoxAlignment -> {"Columns" -> {{Left}}}, 
                  DefaultBaseStyle -> "Column", 
                  GridBoxItemSize -> {
                   "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], 
                 "Column"], "SummaryItem"]}}, 
             GridBoxAlignment -> {
              "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> 
             False, GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
              BaseStyle -> {
              ShowStringCharacters -> False, NumberMarks -> False, 
               PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, 
          GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, 
          GridBoxItemSize -> {
           "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
          BaselinePosition -> {1, 1}]}, 
       Dynamic[Typeset`open$$], ImageSize -> Automatic]},
     "SummaryPanel"],
    DynamicModuleValues:>{}], "]"}],
  SparseArray[
  Automatic, {5, 5}, 0, {
   1, {{0, 2, 5, 8, 11, 13}, {{2}, {1}, {1}, {2}, {3}, {4}, {3}, {2}, {3}, {
     5}, {4}, {4}, {5}}}, {1, -2, 1, -2, 1, 1, -2, 1, 1, 1, -2, 1, -2}}],
  Editable->False,
  SelectWithContents->True,
  Selectable->False]], "Output",
 CellChangeTimes->{3.7850378381362057`*^9},
 CellLabel->
  "Out[103]=",ExpressionUUID->"a88312bd-3924-4d43-8c5b-323efc7aac0e"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"c", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"Dimensions", "[", "A", "]"}], "//", "Flatten"}], "//", 
    "First"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"d", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"A", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "i"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "c"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"u", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"A", "[", 
      RowBox[{"[", 
       RowBox[{
        RowBox[{"i", "-", "1"}], ",", "i"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "2", ",", "c"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"A", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", 
        RowBox[{"i", "+", "1"}]}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", 
       RowBox[{"c", "-", "1"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"i", "=", "2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"While", "[", 
   RowBox[{
    RowBox[{"i", "<", 
     RowBox[{"c", "+", "1"}]}], ",", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"l", "[", 
         RowBox[{"[", 
          RowBox[{"i", "-", "1"}], "]"}], "]"}], "=", 
        RowBox[{
         RowBox[{"l", "[", 
          RowBox[{"[", 
           RowBox[{"i", "-", "1"}], "]"}], "]"}], "/", 
         RowBox[{"d", "[", 
          RowBox[{"[", 
           RowBox[{"i", "-", "1"}], "]"}], "]"}]}]}], ";", 
       RowBox[{
        RowBox[{"d", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "=", 
        RowBox[{
         RowBox[{"d", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "-", 
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"[", 
            RowBox[{"i", "-", "1"}], "]"}], "]"}], "*", 
          RowBox[{"l", "[", 
           RowBox[{"[", 
            RowBox[{"i", "-", "1"}], "]"}], "]"}]}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"b", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "=", 
        RowBox[{
         RowBox[{"b", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "-", 
         RowBox[{
          RowBox[{"b", "[", 
           RowBox[{"[", 
            RowBox[{"i", "-", "1"}], "]"}], "]"}], "*", 
          RowBox[{"l", "[", 
           RowBox[{"[", 
            RowBox[{"i", "-", "1"}], "]"}], "]"}]}]}]}]}], "}"}], ";", 
     RowBox[{"i", "++"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", "=", 
   RowBox[{"c", "-", "1"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"x", "[", 
    RowBox[{"[", "c", "]"}], "]"}], "=", 
   RowBox[{
    RowBox[{"b", "[", 
     RowBox[{"[", "c", "]"}], "]"}], "/", 
    RowBox[{"d", "[", 
     RowBox[{"[", "c", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"While", "[", 
   RowBox[{
    RowBox[{"n", ">", "0"}], ",", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "[", 
        RowBox[{"[", "n", "]"}], "]"}], "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"b", "[", 
           RowBox[{"[", "n", "]"}], "]"}], "-", 
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "n", "]"}], "]"}], "*", 
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"n", "+", "1"}], "]"}], "]"}]}]}], ")"}], "/", 
        RowBox[{"d", "[", 
         RowBox[{"[", "n", "]"}], "]"}]}]}], "}"}], ";", 
     RowBox[{"n", "--"}]}]}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.785037846560337*^9, 3.785037848833924*^9}},
 CellLabel->
  "In[138]:=",ExpressionUUID->"977d2b80-5f41-4207-9ca0-110548410cf6"],

Cell[CellGroupData[{

Cell[BoxData["x"], "Input",
 CellChangeTimes->{3.7850380099755197`*^9},
 CellLabel->
  "In[147]:=",ExpressionUUID->"690cb3dd-81d1-443f-93eb-07e61281d0a7"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox["549", "10"]}], ",", 
   RowBox[{"-", 
    FractionBox["544", "5"]}], ",", 
   RowBox[{"-", 
    FractionBox["1597", "10"]}], ",", 
   RowBox[{"-", 
    FractionBox["2594", "15"]}], ",", 
   RowBox[{"-", 
    FractionBox["757", "6"]}]}], "}"}]], "Output",
 CellChangeTimes->{3.785038010801832*^9, 3.7850380513356566`*^9},
 CellLabel->
  "Out[147]=",ExpressionUUID->"a6a8563f-9b22-47c9-98fb-ae93ecb8525a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["x"], "Input",
 CellChangeTimes->{3.7850379300057993`*^9},
 CellLabel->
  "In[125]:=",ExpressionUUID->"29514d17-6c38-45fa-993a-359b5ca79053"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox["131", "6"]}], ",", 
   RowBox[{"-", 
    FractionBox["128", "3"]}], ",", 
   RowBox[{"-", 
    FractionBox["123", "2"]}], ",", 
   RowBox[{"-", 
    FractionBox["139", "3"]}], ",", 
   RowBox[{"-", 
    FractionBox["157", "6"]}]}], "}"}]], "Output",
 CellChangeTimes->{3.7850379307952356`*^9},
 CellLabel->
  "Out[125]=",ExpressionUUID->"d3d90bdc-9e80-47ac-a1c7-96787f0183b3"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"b", "=", 
  RowBox[{"{", 
   RowBox[{"1", ",", "2", ",", "34", ",", "5", ",", "6"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.785037887905958*^9, 3.7850379042904105`*^9}},
 CellLabel->
  "In[114]:=",ExpressionUUID->"64fb1605-f544-474c-b6ee-f9f939b059ec"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"x", "=", "b"}]], "Input",
 CellChangeTimes->{{3.7850379164131355`*^9, 3.785037918593577*^9}},
 CellLabel->
  "In[115]:=",ExpressionUUID->"7eb212f3-5ef2-4b9f-bc1b-319e91874c29"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "2", ",", "34", ",", "5", ",", "6"}], "}"}]], "Output",
 CellChangeTimes->{{3.785037905721496*^9, 3.7850379192742195`*^9}},
 CellLabel->
  "Out[115]=",ExpressionUUID->"88311c76-d591-4ac6-9f5a-5cfec1b7c804"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Clear", "[", "]"}]], "Input",
 CellChangeTimes->{{3.7850378544124885`*^9, 3.785037856563903*^9}, {
  3.7850378992342944`*^9, 
  3.7850379001564198`*^9}},ExpressionUUID->"d91148fd-1f6c-4042-8dc9-\
bf0a2f17d9ad"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.7850389758133864`*^9, 
  3.785039050567207*^9},ExpressionUUID->"438ae34c-3aba-41ef-9f22-\
517d9c4df33c"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1280, 637},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
CellContext->Notebook,
FrontEndVersion->"11.3 for Microsoft Windows (64-bit) (2018\:5e743\:670828\
\:65e5)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 186, 3, 96, "Title",ExpressionUUID->"8961cf75-fdc5-4109-b21b-7a5dddeaa7ee"],
Cell[769, 27, 215, 3, 45, "Input",ExpressionUUID->"499f085a-2c33-40be-ba39-97a7fef4f643"],
Cell[987, 32, 8686, 226, 867, "Input",ExpressionUUID->"39468ca8-ce8d-4ece-893c-13e29e110352"],
Cell[CellGroupData[{
Cell[9698, 262, 240, 4, 45, "Input",ExpressionUUID->"2139ec7d-9357-4f15-8c94-601722cbade6"],
Cell[9941, 268, 922, 29, 103, "Output",ExpressionUUID->"b1e4a0d4-affe-4465-a7ab-1fba7e15a44e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10900, 302, 246, 4, 45, "Input",ExpressionUUID->"125fc4cd-f5fb-4769-b181-a7e7a26ef470"],
Cell[11149, 308, 922, 28, 103, "Output",ExpressionUUID->"82492a5e-7ab0-49a8-934a-53fbd03ef124"]
}, Open  ]],
Cell[12086, 339, 11522, 266, 848, "Input",ExpressionUUID->"03cf18a0-f1e3-478c-af1d-32b85ea495ba"],
Cell[23611, 607, 483, 12, 28, "Input",ExpressionUUID->"9a0bdda5-fff6-4233-8d25-85399e077918"],
Cell[CellGroupData[{
Cell[24119, 623, 204, 3, 45, "Input",ExpressionUUID->"ab511fbc-18fa-4bdc-acea-d91f643b0124"],
Cell[24326, 628, 757, 22, 79, "Output",ExpressionUUID->"b26c7702-0cba-40fb-b94c-72b1b8c3a1a7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25120, 655, 150, 2, 28, "Input",ExpressionUUID->"c2fc334d-9952-414a-ba35-3b4dfd2b27f9"],
Cell[25273, 659, 476, 11, 32, "Output",ExpressionUUID->"7b46d1ba-0090-4ece-8357-9e57b7b2b8a5"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25786, 675, 297, 5, 45, "Input",ExpressionUUID->"51d13254-0df7-4704-bf49-d6286f85ad74"],
Cell[26086, 682, 1082, 27, 79, "Output",ExpressionUUID->"a0f1a3bf-4e2d-41f0-9bba-91e7304a71bd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27205, 714, 206, 3, 45, "Input",ExpressionUUID->"663a59d5-690f-4c7d-8e78-d9e8a1a18fec"],
Cell[27414, 719, 757, 22, 79, "Output",ExpressionUUID->"897ba509-0957-4b4d-aabe-113893b5f538"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28208, 746, 670, 18, 45, "Input",ExpressionUUID->"734c34da-98a4-4a5a-9511-b66e4e4a216d"],
Cell[28881, 766, 173, 2, 32, "Output",ExpressionUUID->"12bb7bd1-88b0-469b-8c2f-679f0613800a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29091, 773, 280, 8, 45, "Input",ExpressionUUID->"2b64211b-bf82-4c8f-863d-a458f87781e5"],
Cell[29374, 783, 498, 14, 32, "Output",ExpressionUUID->"98129c92-b617-497c-8e55-169e764a3ed6"]
}, Open  ]],
Cell[29887, 800, 219, 5, 45, "Input",ExpressionUUID->"187d9754-897c-4f34-84ae-95937a6d96d5"],
Cell[30109, 807, 313, 8, 28, "Input",ExpressionUUID->"2513bbad-9528-41ae-b06c-8ca0a2ec8c7b"],
Cell[30425, 817, 206, 4, 32, "Output",ExpressionUUID->"f520807b-dc45-476b-bbbb-79d9c13c130c"],
Cell[CellGroupData[{
Cell[30656, 825, 242, 4, 45, "Input",ExpressionUUID->"8ebc8e24-1a22-44f6-a9ad-b90763b20b01"],
Cell[30901, 831, 1576, 49, 144, "Output",ExpressionUUID->"e7a72ae7-a23c-40a2-b4e1-47cdee59c48c"]
}, Open  ]],
Cell[32492, 883, 329, 10, 8, "Text",ExpressionUUID->"c3562191-a342-44d0-842e-b936de6643f4"],
Cell[CellGroupData[{
Cell[32846, 897, 240, 4, 67, "Section",ExpressionUUID->"544df9ba-2d11-4896-b912-9d4bb010df07"],
Cell[33089, 903, 8932, 231, 737, "Input",ExpressionUUID->"3d281545-cce5-40a1-9b8c-9d3906fa76e4"],
Cell[CellGroupData[{
Cell[42046, 1138, 150, 2, 28, "Input",ExpressionUUID->"eb3ac092-d10c-40cf-ae97-ef6901f6f0ff"],
Cell[42199, 1142, 499, 12, 32, "Output",ExpressionUUID->"9d475d11-c9cd-41dc-a54d-8ab9b738409a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42735, 1159, 431, 10, 45, "Input",ExpressionUUID->"5d92213f-562b-4370-b4b9-c86987497bd4"],
Cell[43169, 1171, 1058, 30, 103, "Output",ExpressionUUID->"116b9dae-6cad-4a0d-8ce0-62b664e6f91f"]
}, Open  ]],
Cell[CellGroupData[{
Cell[44264, 1206, 147, 2, 28, "Input",ExpressionUUID->"a953f32c-a9d1-4965-b2e4-be2bea9a1512"],
Cell[44414, 1210, 381, 10, 32, "Output",ExpressionUUID->"0f144c48-9519-4e5c-9345-213c7cf6ef57"]
}, Open  ]],
Cell[44810, 1223, 5159, 138, 531, "Input",ExpressionUUID->"44bdc003-a7c7-4e48-ba61-958730091d5a"],
Cell[49972, 1363, 329, 10, 8, "Text",ExpressionUUID->"931b6f6b-47e3-4ed2-b7f2-f6c41a4667e7"],
Cell[50304, 1375, 329, 10, 8, "Text",ExpressionUUID->"6db80a03-209f-42e1-9d49-8435f4a92963"]
}, Open  ]],
Cell[CellGroupData[{
Cell[50670, 1390, 165, 3, 67, "Section",ExpressionUUID->"ab435dfa-bde9-4bf5-975d-3e57cf37e8da"],
Cell[CellGroupData[{
Cell[50860, 1397, 265, 5, 53, "Subsection",ExpressionUUID->"0398b880-538e-4f2e-8967-143e7eec2da6"],
Cell[51128, 1404, 1337, 36, 119, "Input",ExpressionUUID->"01c3cd34-de03-42fa-a4a2-72ef6acc146b"],
Cell[52468, 1442, 329, 10, 8, "Text",ExpressionUUID->"4918d143-afdc-472a-924b-c798cf3f2ccd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52834, 1457, 276, 5, 53, "Subsection",ExpressionUUID->"10b3b4bc-2a60-419f-8c69-026bc319969d"],
Cell[53113, 1464, 3569, 84, 290, "Input",ExpressionUUID->"99913f9a-9c8f-4532-8019-26c7c1448f3a"],
Cell[56685, 1550, 329, 10, 8, "Text",ExpressionUUID->"3ae9cd2e-261e-4529-9615-c6ae4663ea64"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[57063, 1566, 165, 3, 67, "Section",ExpressionUUID->"1c2529f6-eba0-49be-a041-cf15b9ec240f"],
Cell[57231, 1571, 5310, 143, 436, "Input",ExpressionUUID->"9c2dd963-5d83-495a-8b25-173af639ecfa"],
Cell[CellGroupData[{
Cell[62566, 1718, 642, 20, 45, "Input",ExpressionUUID->"03fcadc4-b839-4aeb-be71-dd81d9e7a5b5"],
Cell[63211, 1740, 9083, 191, 71, "Output",ExpressionUUID->"1cbcbc5e-96c0-4c32-a8e5-63a0e47971ad"]
}, Open  ]],
Cell[72309, 1934, 272, 6, 45, "Input",ExpressionUUID->"3d552d3b-61f1-487d-b536-b3c8bf549289"],
Cell[CellGroupData[{
Cell[72606, 1944, 156, 3, 28, "Input",ExpressionUUID->"20d58777-a045-42db-9dcb-e7d954b5f39b"],
Cell[72765, 1949, 151, 2, 32, "Output",ExpressionUUID->"b6d4b11d-a8ea-45a1-abd0-aa4be3eaeb95"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72953, 1956, 275, 7, 28, "Input",ExpressionUUID->"81b5ab71-30a3-4ea0-9925-12ca72280e40"],
Cell[73231, 1965, 179, 3, 32, "Output",ExpressionUUID->"abb9fca9-028a-4a19-8e0b-8912d0fd8521"]
}, Open  ]],
Cell[CellGroupData[{
Cell[73447, 1973, 318, 8, 28, "Input",ExpressionUUID->"b2987884-cf18-49fd-8c43-a84c7dc5f765"],
Cell[73768, 1983, 476, 15, 55, "Output",ExpressionUUID->"61608589-1733-4257-b2e0-468b406b31c7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[74281, 2003, 317, 8, 28, "Input",ExpressionUUID->"5a6c3f87-f334-4d9f-be45-aa6e1747f48b"],
Cell[74601, 2013, 847, 21, 55, "Output",ExpressionUUID->"3775de8f-cad6-4103-8547-2041f887d7bb"]
}, Open  ]],
Cell[CellGroupData[{
Cell[75485, 2039, 200, 4, 28, "Input",ExpressionUUID->"5df2f2be-d24c-4e19-ab6b-b64e1d9914e9"],
Cell[75688, 2045, 9089, 192, 71, "Output",ExpressionUUID->"a88312bd-3924-4d43-8c5b-323efc7aac0e"]
}, Open  ]],
Cell[84792, 2240, 3853, 126, 327, "Input",ExpressionUUID->"977d2b80-5f41-4207-9ca0-110548410cf6"],
Cell[CellGroupData[{
Cell[88670, 2370, 154, 3, 28, "Input",ExpressionUUID->"690cb3dd-81d1-443f-93eb-07e61281d0a7"],
Cell[88827, 2375, 477, 15, 55, "Output",ExpressionUUID->"a6a8563f-9b22-47c9-98fb-ae93ecb8525a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[89341, 2395, 154, 3, 28, "Input",ExpressionUUID->"29514d17-6c38-45fa-993a-359b5ca79053"],
Cell[89498, 2400, 450, 15, 55, "Output",ExpressionUUID->"d3d90bdc-9e80-47ac-a1c7-96787f0183b3"]
}, Open  ]],
Cell[89963, 2418, 278, 6, 28, "Input",ExpressionUUID->"64fb1605-f544-474c-b6ee-f9f939b059ec"],
Cell[CellGroupData[{
Cell[90266, 2428, 200, 4, 28, "Input",ExpressionUUID->"7eb212f3-5ef2-4b9f-bc1b-319e91874c29"],
Cell[90469, 2434, 255, 5, 32, "Output",ExpressionUUID->"88311c76-d591-4ac6-9f5a-5cfec1b7c804"]
}, Open  ]],
Cell[90739, 2442, 234, 5, 45, "Input",ExpressionUUID->"d91148fd-1f6c-4042-8dc9-bf0a2f17d9ad"],
Cell[90976, 2449, 152, 3, 28, "Input",ExpressionUUID->"438ae34c-3aba-41ef-9f22-517d9c4df33c"]
}, Open  ]]
}, Open  ]]
}
]
*)

